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The Pauli exclusion principle is a constraint on the natural occupation numbers of fermionic states. It has been 
suspected since at least the 1970's, and only proved very recently, that there is a multitude of further constraints 
on these numbers, generalizing the Pauli principle. Here, we provide the first analytic analysis of the physical 
relevance of these constraints. We compute the natural occupation numbers for the ground states of a family of 
interacting fermions in a harmonic potential. Intriguingly, we find that the occupation numbers are almost, but 
not exactly, pinned to the boundary of the allowed region (quasi-pinned). The result suggests that the physics 
behind the phenomenon is richer than previously appreciated. In particular, it shows that for some models, the 
generalized Pauli constraints play a role for the ground state, even though they do not limit the ground-state 
energy. Our findings suggest a generalization of the Hartree-Fock approximation. 
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Introduction. — In 1925 the study of atomic transitions led 
to Fault's exclusion principle |1|. It states that for identical 
fermions the occupation number for any quantum state cannot 
exceed the value 1. By 1926, Dirac [2 J and Heisenberg |3| had 
identified the exclusion principle as a consequence of a much 
deeper statement: the anti-symmetry of the many-fermion 
wave function. While anti-symmetry allows one to find the 
correct solutions to the full many-fermion Schrodinger equa- 
tion, it does not render the exclusion principle obsolete: In 
many situations, the latter is sufficient to predict the qualita- 
tive behavior of fermionic systems without the need to resort 
to (often computationally intractable) ab initio methods. The 
Aufbau principle for elements in the periodic table serves as a 
prime example. 

This observation motivates the study of generalizations of 
the exclusion principle, which, maybe surprisingly, exist and 
exhibit an extremely rich structure |4). To set the scene, note 
that the Pauli constraint can be stated succinctly as 

< Ai < 1 Vi (1) 

in terms of the natural occupation numbers Xi, which are 
the eigenvalues of the 1-particle reduced density operator 
(1— RDO) pi, normalized to the particle number N . The util- 
ity of the exclusion principle is grounded in the fact that in the 
ground states of many-fermion systems, one often observes 
Ai « or A; « 1, which is equivalent to stating that the 
Hartree-Fock approximation works fairly well in these sys- 
tems. 

It had been observed in the 1970s that there are further lin- 
ear inequalities respected by the natural occupation numbers 
as a result of global anti-symmetry ISIQ. One particular ex- 
ample is the so-called Borland-Dennis setting A"^ [H,q\ of three 
fermions and a six dimensional 1— particle Hilbert space T-Lq 
||2J. Here, the set of constraints is given by 

Ai + Ae - A2 + Ag = A3 + A4 = 1 (2) 
:= A5 + Ae - A4 > (3) 

on the ordered eigenvalues Ai > A^+i. 



In a ground-breaking work building on recent progress in 
invariant theory and representation theory, Klyachko exhib- 
ited an algorithm for computing all such Pauli-like constraints 
pl, f8l|. In fact, his work is part of a more general effort in 
quantum information theory addressing the quantum marginal 
problem which asks when a given set of single-site reduced 
density operators (marginals) is compatible in the sense that 
they arise from a common pure global state (see also 19]- 
[T2I ). The global state may be subject to certain symme- 
try constraints — one obtains the fermionic case (commonly 
known as the iV— representability problem iTOlfl?! ') by requir- 
ing total anti-symmetry under particle exchange. Klyachko 
showed that for fixed particle number N and dimension d of 
the 1— particle Hilbert space, the generalized Pauli constraints 
amount to affine inequalities of the form 

KO + KlAl + . . . + K^Xd > 0. (4) 

Geometrically, these constraints define a convex polytope 
'PN,d C M'' of possible spectra (for more details see Appendix 
[A| . In general, if a spectral inequality such as ([T]l or Q is 
(approximately) saturated we say that the corresponding spec- 
trum is (quasi-) pinned to its extremum. 

The natural question arises as to whether ground states of 
relevant many-body models saturate some of those inequal- 
ities. Strong numerical evidence supporting this conjecture 
has been presented in (15). The problem is challenging to ad- 
dress analytically, as one has to not only compute the ground 
state, but also determine and diagonalize the corresponding 
1 RDO. 

Here, we present for the first time an analytic analysis. For 
the ground state of a model of interacting fermions in a har- 
monic potential, the natural occupation numbers are calcu- 
lated. We obtain several results. We confirm that for this 
very natural model, the natural occupation numbers lie indeed 
close to the boundary of set of allowed ones. The analytic 
analysis enables us to track the "trajectory" of eigenvalues as 
a function of the interaction strength between the fermions. 
What is conceptually also important, is the fact that the eigen- 
values never lie exactly on the boundary. To see why one 
could expect the opposite, note that the ground state energy of 
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a Hamiltonian H = ^ /i^*^^) with two-particle terms /i^*'^) 
can be represented as a constrained optimization problem 



N 

E^in = mill V tr[/i(* 



where the P2 are 2-particle density operators that are com- 
patible in the sense that they are the reduced densities of some 
A^— fermion state [13, 14). Since the energy functional is lin- 
ear, it does not possess an unconstrained minimum. There- 
fore, E'niin will be achieved on the boundary of the set of com- 
patible density operators, where at least one of the compati- 
bility constraints is active in the sense that any further min- 
imization would violate it. One way of understanding why 
a "pinning" effect for the natural occupation numbers is ob- 
served, is to posit that the generalized Pauli constraints are 
among the active physical constraints. While this effect may 
well occur, we show in this work that q'Mas/-pinning appears 
in natural fermionic systems: the eigenvalue constraints seem 
to play a role, but are not active in the above sense. The find- 
ing suggests that the physics of the phenomenon is richer than 
previously appreciated. We will return to the physical conse- 
quences quasi-pinning has on the structure of wave functions 
after presenting the calculations for our model system. 

The Model. — In order to analyze possible pinning effects 
analytically, we consider a model of N identical fermions sub- 
ject to a harmonic external potential and a harmonic interac- 
tion term: 
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(5) 



The corresponding eigenvalue problem without any symmetry 
constraint can easily be solved by transforming the Hamilto- 
nian to the one of decoupled harmoiuc oscillators. Two eigen- 
frequencies appear: a non-degenerate one describing the cen- 
ter of mass motion and another {N — l)-fold degenerate fre- 
quency associated with the relative motion. The natural length 
scales corresponding to these eigenmodes are 
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muj ' 



mcj^l + ND/{muj'^)' 



By rescaling the energy and the length scale, the fermion- 
fermion coupling constant D can be absorbed by the term 
moj^. Hence the spectrum A of a 1— RDO corresponding to an 
eigenstate of H depends only on the relative fermion-fermion 

interaction strength ~ (l) ~ 1- I" fact, it will prove 
slightly more convenient to parameterize the coupling using 



6 := In 



04- 



ND 



(6) 



Then, in the regime of weak interaction, D and 6 are in leading 



order proportional, D = 



N 



To study the physical relevance of the generalized Pauli 
constraints we restrict the Hamiltonian H to the fermionic 



Hilbert space A^['Hoo], with Tioo = i (K), i-e. we are treat- 
ing the N particles as fermions (without spin). In fW\ H has 
been diagonalized and the ground state reads in spatial repre- 
sentation {x — (xi, . . . , Xn)) 



vI'jy(x) = const X 

l<i<j<Ar 

2N \P p 



Xj) 



(7) 



X exp 
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(Note its structural similarity to Laughlin's ground state wave 
function describing the fractional quantum Hall effect |17|. 
Moreover, the polynomial in front of the exponential function 
is the Vandermonde determinant and by omitting it we obtain 
the ground state in the bosonic iV-particle Hilbert space.) 

The spectrum and its properties. — We now outline the 
calculation of the spectrum X{5) as a function of the coupling. 
We omit details of this tedious but mostly straight-forward 
computation, presenting the final result, together with some 
conceptual insights obtained along the way. 

The 1— RDO is calculated by integrating out — 1 coor- 
dinates of the A^-fermion state pn{x,x') = \E'^(x)\E'Ar(x'). 
An exercise in Gaussian integration and integration by parts 
yields 

Pi{x, x) = p{x, x) exp[— a(x^ + x"^) + j3xx] , 



where p is a symmetric polynomial of degree (^) in the vari- 
ables X, x' originating from the Vandermonde determinant in 
(|7]i, and a and /3 some constants depending oxvl,! and N . 

If the fermions do not interact with each other, the ground 
state l^'jv) is a single Slater determinant and the spectrum of 
its 1— RDO is trivial, i.e. 



A((5 = 0) - (1,...,1,0,...). 



(8) 



JV 



The regime of weak interaction can be characterized by the 
condition |D| ^ rrnxP' or, equivalently, (5 w 0. We thus 
employ degenerate perturbation theory to obtain A((5) around 
(5 = 0. The reason we employ the parameter b is that one can 
prove a duality 



Afe (5) = Afc {-b) 



(9) 



relating the spectra for attractive (8 < 0) and repulsive {S > 0) 
fermion-fermion interaction (interestingly, that this duality 
holds is not obvious on the level of ground-state wave func- 
tions). This immediately implies that the expansion X{S) con- 
tains only even order terms, simplifying the perturbation the- 
ory. 
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The solution for = 3 reads: 
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Similar results follow for N 
chy of the eigenvalues. 



(10) 

2. Note the non-trivial hierar- 



(11) 



for all fc > 5. Moreover, the spectrum A for 5 not too large 
is very close to the one of a single Slater determinant. For 
instance, A^, i = 1,2,3 deviate from 1 and Aj, j > 4 from 
only by at most 1 percent if \6\ < 0.5. This emphasizes the 
relevance of the Pauli constraints ([T]i. 

Quasi-Pinning by Generalized Pauli Constraints. — 
Equipped with the explicit solution ( [T0| , we can proceed 
to analyze whether the generalized Pauli constraints play a 
role for the ground state. While the underlying 1— particle 
Hilbert space Hoo is infinite-dimensional, the scaling (111 
implies that the spectrum is strongly concentrated on a 
low-dimensional subspace, at least for small S. One can use 
this fact to deduce statements about the position of the total 
eigenvalues from truncated information alone. 

This can be understood from simple geometric consid- 
erations. Let d < d' < oo. Because a d-dimensional 
1— particle Hilbert space can be imbedded into any (larger) 
d'-dimensional one, one sees that the convex polytope 'PN,d is 
nothing but the intersection between Vna' ^nd the set of spec- 
tra with only d non-zero eigenvalues (see also Appendix [BJ. 
Hence any facet of Vn.iI arises from the intersection of some 
facet of VN.d' with the subspace of said spectra. Formally, a 
facet F' of Vn.j.' consists of points saturating a generalized 
Pauli constraint 



2 = 1 



lAi + 



d' 

E ' 



,X > 0. 



(12) 



Denote the first two summands by Z?(A*''), where A*'' ~ 
{^i)i=i is the truncated spectrum. Clearly, £'(A*'') — de- 
scribes the restriction of the facet to the d-dimensional setting. 
Now assume the truncated spectrum A*'' ((5) is not pinned, i.e. 
£'(A*' ((5)) > 0, then the hierarchical scaling i 
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implies 



2d-4\ 



which is positive for 5 small enough. Hence the full spectrum 
A' also fails to be pinned. The case d' — oo works in the same 
way, up to some mild assumptions (see Appendix [BJ. 

We will now apply these considerations to our model. First, 
we truncate to 6 dimensions, which has the advantage that the 
spectral polytope corresponding to A'^['H6] is 3-dimensional 
and can thus be visualized. In a second step, we take a seventh 
eigenvalue into account. This setting turns out to be strong 
enough to establish all statements we have mentioned above - 
namely that the total spectrum is not exactly pinned, but does 
lie close to the boundary (quasi-pinned). 
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(o) ■ 



FIG. 1: Spectral "trajectory" v{S) (thick line, partially covered by 
facet, schematic) up to correction of order and small part of 
the polytope V around vertex u'"' obtained by cutting V along the 
dashed lines 

The simplest non-trivial setting /\^[He] becomes an appro- 
priate description if A7, Ag, . . . ~ 0. By ( [TT| , this condition 
is fulfilled if 6 is small enough that contributions of order 
can be neglected. Choosing A4, A5 and Ag as free parameters 
according to (|2ji, the corresponding polytope 7^3 g effectively 
reduces L15J to a 3-dimensional polytope V C M.^ with ver- 
tices. 



w*"^ = (0,0,0) 

-/ill 
2' 4' 4 



(b) 



' 2' 2' 2 



(14) 



Hence the vertex v^"-' corresponds to single Slater determi- 
nants and the 2-facet spanned by w^"), w^'') and w*^"^^ is defined 
by 13*^^^ = 0, which is here the one of interest and represents 
exact pinning by constrain t ((3] l. We first illustrate schemati- 
cally our result ( [T0| in Fig.[T] There, the spectral "trajectory", 

is shown as a thick line (neglecting effects of order 5^ and 
higher). It starts at the vertex w'"^ which corresponds to the 
non-interacting situation (5 = 0. When increasing the fermion- 
fermion interaction, v{S) leaves the vertex and moves 
along the edge (u^''^ u^''-'), the distance to v^"-'> growing as 
6^. On the finer scale S^, v{6) also moves away from the edge 
but is still pinned to the boundary of the polytope, lying on 
the 2-facet spanned by u^"', i;^'') and v'^'^\ This is the bottom 
area in Fig.[T] corresponding to the constraint ([3]). 

The pinning seems to disappear if we consider higher or- 
ders. From ( [TOj i, we can infer that the distance to the 2-facet 
(y(o-) ^ tiC*) J increases as 5*, 



(13) 



(15) 
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with C*^^^ — ^^jjf). However, this calculation is inconclusive, 
as the distance to the boundary is of the same order, (5®, as the 
truncation error (recall (fTSjl). 

To resolve the issue, we take another eigenvalue, A7, into 
account. We thus work in the setting A'^['H7] with four con- 
straints D^p > for i = 1, . . . , 4 [4|. This setting is valid as 
long as As , Ag , . . . « or in other words we neglect terms of 
order (5^° or higher (but in contrast to the setting A'^['H6] we 
include 5*— terms). Since the polytope is now 6— dimensional 
we cannot present our results graphically anymore. The re- 
sults ([10} lead to a = 1,2,3,4) 

dP = CPS^ + 0(<5iO) , (16) 

witVi ^C^) - r^^) - 20_ A'^') — ^<L A"^^ — 2890 

Willi — 2187, 1,2 — 243' ^3 ~ 2187' ''4 — 59049' 

Here in the [■H7] -analysis, the new result is that all four dis- 
tances Dp' are non-zero to a smaller order, 6^, than the error 
of spectral truncation, (5^°. Together with the comments at 
the beginning of this section, this shows that the absence of 
pinned spectra is genuine, rather than an artifact of the trunca- 
tion. Given this, the quasi-pinning found here is surprisingly 
strong. In particular it exceeds by four additional orders the 
(quasi-)pinning by Pauli's exclusion principle constraints ([T]l, 

< 1-A2(5),1-A3((5),A4((5),A5(<5) = ^SU0{6^). (17) 

Generalizing Hartree-Fock. — In this section, we dis- 
cuss what conclusions can be drawn about the iV— fermion 
state \^>) itself, given information just about the position 
of the eigenvalues of the corresponding 1-RDO relative to 
the boundary of the spectral polytope. In this way, quasi- 
pinned spectra are endowed with a physical significance. 
To this end, recall the basic fact that the spectrum Agi = 
(1, . . . , 1, 0, . . . , 0) can only arise from a Slater determinant 
\'^) — |1, . . . , TV). It is well-known that this statement is 
stable under small deviations: if A sa Asi, then I^I^) is well- 
approximated by a Slater determinant (see [18] or Appendix 
[D). 

For exactly pinned spectra, there is a simple generalization 
of these observations. In ifTSl . it is stated that constraint (|3]l 
can be saturated only by states of the form 

I*) =a|l,2,3)+^|l,4,5)+7|2,4,6), 

a fact is dubbed "selection rule" for Slater determinants (see 
also Appendix [c]). The general statement reads: if D{X) > 
is a generalized Pauli constraint, then D{\) — can only 
be achieved by states jvP) which are superpositions of those 
Slater determinants whose (unordered) spectra also saturate 
D. 

What is more important, a stable version of this state- 
ment applying to quasi-pinned states can be found — at least 
for specific situations. In the Appendix [D] we show that 
for the Borland-Dennis setting, spectra in the vicinity of the 
facet corresponding to constraint ^ are approximately of the 
form above. In particular, quasi-pinned states are close to 
states containing fairly low amounts of multi-partite entan- 
glement as quantified by the Schmidt number U9J . In EOl 



a new entanglement measure has been suggested, which, for 
the Borland-Dennis setting, naturally separates exactly pinned 
and non-pinned states. We believe that these findings open 
up a potentially significant avenue for investigating the struc- 
ture of fermionic ground states via their natural occupation 
numbers — generalizing a program that has long been carried 
out for the Hartree-Fock case 1 1 8 1 . 

We close by speculating that these insights could give rise 
to improved numerical procedures. The idea is to replace the 
ground state ansatz of one single Slater determinant by the 
states corresponding to the points lying on the (quasi-)pinning 
polytope facet. In contrast to the configuration interaction (CI) 
methods in quantum chemistry which improve the Hartree- 
Fock approximation by adding several arbitrary Slater deter- 
minants to the Hartree-Fock state our method would add only 
a few but carefully chosen additional Slater determinants. 

Conclusions. — For a natural model of interacting 
fermions in a harmonic trap we analytically calculated the 
leading orders of the eigenvalues of the 1— RDO correspond- 
ing to the fermionic ground state as function of 5, a measure 
for the fermion-fermion interaction strength. The investiga- 
tion of the generalized Pauli constraints has shown that none 
of them is completely saturated, which might be a generic 
property of all continuous models of interacting fermions. In 
particular, the findings show that it is likely extremely chal- 
lenging to use numerical methods to distinguish between gen- 
uinely pinned and mere quasi-pinned states. This underscores 
the need for analytical analyses, first provided here. On the 
other hand the pinning up to corrections of order 6^ we found 
here is surprisingly strong. In particular it exceed the one by 
the Pauli exclusion principle constraints ([T]i, which are pin- 
ning up to corrections of order S"^ only. 
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Appendix 

This appendix is split into four sections. The first one in- 
troduces the notation and repeats the solution of the fermionic 
quantum marginal problem. In the second section we explain 
how to simplify the pinning analysis by truncating the spec- 
trum. This amounts to the proof of statement (13), a relation 
connecting polytope distances of the correct and truncated 
marginal setting. The third section introduces a selection rule, 
which explains how the structure of a A^— fermion state sim- 
plifies if its natural occupation numbers are exactly pinned to 
some Pauli facet and applies it to the Borland-Dennis setting. 
In the last section we present a modification of this selection 
rule for the case of only approximate pinning. This then justi- 
fies our Hartree-Fock generalization. 
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A. Notation and Fermionic Quantum Marginal Problem. — 



The problem of determining all spectra 



A = (A.)ti 



1 > Ai > A2 > . . . > Ad' > 



(18) 



of 1— particle reduced density operators (1— RDO) pi arising 
from some pure A^— fermion state \^) G A^['Hd'], 



(19) 



by tracing out — 1 particles, is known as the fermionic quan- 
tum marginal problem of the setting l\^[Hd']- Here d' E 
N U {00}, Hd' is the d'— dimensional separable 1— particle 
Hilbert space and we use the trace normalization convention, 



ti-[pi] = Ai 



.. + Xd'=N, 



(20) 



common in quantum chemistry. 

For d' finite, the family of possible spectra (we call them 
compatible w.r.t A^['Hd']), is described by finitely many in- 
dependent conditions {C;}, the generalized Pauli constraints. 
Each of them has the form 



Cj : A(A) = Ko + KiXi + ... Kd'Xd' > 0, 



(21) 



kq, . . . ,Kd' G Z and describes a half-space ofW^'. These 
constraints together with the trivial conditions ( fTsj l and ( 20 1 
define the polytope VN.d' C M'' of possible spectra. In that 



sense every constraint (21 1 gives rise to a facet Fi of this poly- 
tope. 



F, = {XeVN,d'\D^{X)^0}■ 



(22) 



Note that besides these Pauli facets there are also further 
facets, those corresponding to the trivial constraints ( fTSj l, but 
they will not be of interest in our work. Moreover, the quantity 
Di{-), which is only defined up to a positive factors, defines 
after fixing this factor (i.e. the parameters Ki) a measure for 
the distance of spectra to the corresponding facet Fi . For the 
case of d' finite it coincides up to normalization with the Eu- 



clidean distance, dist2(/i, Fi) 



, K 



(Hi, 



, Kd'). 



For the case d' — 00 the set Vn.oo of compatible spectra 
is not explicitly known yet. Nevertheless, for our work we 
assume that it is also defined by a family of linear inequalities 



d\°°\x) = no + KiXi + K2X2 + 



> 0. 



(23) 



The results on truncation of the spectrum and the relation of 
polytope VN.d and VN,d', d < d' finite presented in Appendix 
[B]strongly emphasizes that this assumption is justified. More- 
over, the involved fact that the closure Vn^ is convex also 
suggests this assumption. 

Finally, we still make some comments on the meaning of 
natural orbitals {|fc)}, the eigenvectors of the 1— RDO, 



d' 



= ^ Afe|fc)(fc|, 



(24) 



and their utility for applications. 

These natural orbitals induced by a fixed state G 
A^[-Hd'], rf' e N U {00}, define a basis Bi := 
for the 1— particle Hilbert space Hd'- For ease of notation 
we skip the argument 5* of \i{'i>)). Basis Bi then induces the 
basis Bn for A^[Hd'] of corresponding Slater determinants 
(I < ii < . . ■ < iN < d') 



= A 



n[\h) 



\iN)], (25) 



where An is the anti-symmetrizing operator on the 
A^— particle Hilbert space Hj^ ■ expanding | 'if) w.r.t. to 
Bn, 



the natural occupation numbers are given by 

i, kei 



(26) 



(27) 



To compare marginal settings of different dimensions, d, d' 
with d < d' < 00 we imbed Hd into Hd', 



span{|z)}ti ^Hd< Hd' = span{|i)}ti, 



(28) 



where the closure is only relevant for the case d' infinite. In 
the same way. 



a"" [Hd]<A''[Hd']. 
Indeed, according to ([26]l, we find that every state 



I*) 



E 



Ci 1 



eA^[Hd] 



(29) 



(30) 



l<ii < . ..<iN<d 

can be imbedded into A^[l-Ld'] by 



l<2i < ...<iN<.d 



eA''[Hd' 



(31) 



and all the other coefficients Ci in ([31), those with in > d, 
vanish. We used here different symbols for the states I^I^) and 
1 4*') to distinguish between the two different spaces A^\Hd\ 
and A^[Hd'] to which they belong. This subtle difference 
is becoming relevant if we determine the natural occupation 
numbers A' of |*') (recall (p7|). 



(32) 



A' = (Ai,...,A<j,0,...,0) 

d'-d 



fe=l 



differing from A = (Ai, . . . A^) by additional zeros. In the 
following, to simplify the notation, we will use the same sym- 
bols for mathematical objects and their imbeddings into larger 
spaces. 
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B. Truncation of the Spectrum. — 



In our work we have determined the "trajectory" of spectra 



m = {Ms))r=i e ^3,00, 



(33) 



of the 1— RDO corresponding to the ground state of a 
3— fermion model with relative fermion-fermion interaction 
strength S. The goal was then to show that for S not too 
large, X{6) is almost but not exactly saturating some of the 
generalized Pauli constraints of its setting A^[?^oo]- Geomet- 
rically this means that the vector X{S) is very close to some 
Pauli facet Fi of Vs.oo- In that case we say that the spectrum 
is quasi-pinned to the facet Fi. Since 7^3,00 is not explicitly 
known and quite involved (it is described by infinitely many 
constraints on infinitely many eigenvalues), we have truncated 
the spectrum and simplified the pinning analysis by consider- 
ing only the largest d eigenvalues. 



A*' — (Ai, . . . , Xd), 



(34) 



and analyzed the saturation of the constraints corresponding 
to the setting /\'^[Hd]- The following fact justifies this ap- 
proach: For d < d' every Pauli facet F of VN,d is contained 
in some Pauli facet F' of VM,d', i-e. F is the intersection of 
F' with the hyperplane of spectra with only d non-zero eigen- 
values. Then, for small A^+i, Xd+2^ ■ ■ ■, small distance of A*'' 
to F translates to small distances of A to F' modulo an error 
of order of the largest neglected eigenvalue, Xd+i- To illus- 
trate this, we present the example A'^['H6], which is one of the 
two settings studied in our work. There one generahzed Pauli 
constraint reads iHHH 

Z?(6)(A):=2-(Ai+A2 + A4)>0. (35) 

For the setting [Hqo] the known constraint j?) 

(A) = 2 - (Ai + A2 + A4 + A7 + An + A16 + . . .) > , 

(36) 



coincides with constraint ( 35 1 up to a linear combination of 
eigenvalues A7, An, Aig, • • ■, which where neglected in the 
truncated setting. 

A first important step in proving the universality of this re- 
lation between polytope distances of correct and truncated set- 
ting is the next lemma: 

Lemma 1. Consider the quantum marginal problems of the 
two settings A^[H<i] and A^['Hd'], d < d' € N U {00} and 
let X — (Ai, . . . , Xd) be a spectrum. Then, 

(Ai, . . . , Xd) compatible w.r.t. [Hd] 



(Ai, . . . , Xd, 0, . . . , 0) compatible w.r.t. [Hd' 



d'-d 



For the corresponding polytopes this means 



(37) 



(38) 



the polytope V^.d' intersected with the hyperplane given by 
Xd+i, A(i+2, . . . = coincides with VN.d- 



Proof. The direction was already explained at the end of 
Sectio i(a] To prove "<^=" we show that a state \'^') expanded 
according to ( 26 1, 



l<il<...<2jV<'^' 



Ci|i> eA^[Hd'], 



(39) 



with natural occupation numbers (Ai,...,Aci,0,...,0) con- 



tains only Slater determinants with ii 



this is clear due to ( 27 1, which then yields 



yk > d : = Ai 



E 

i, fcei 



Hence Cj = if i at > d. 



, i < d. But 



(40) 



□ 



What does Lemma [T] imply for the relation between the 
families of generalized Pauli constraints of two settings? Let 
us consider two settings with d, d' finite, d < d'. Every con- 
straint Dj for the setting A^[Hd'] is linear and hence its re- 
striction 

^'(Ai, 



,Xd) = D'JXu...,Xd,0,...)>0 



(41) 



to the hyperplane defined by = A^+i, Xd+2, ■ • . is also a lin- 
ear constraint in the remaining coordinates Aj, . . . , A^. How 
is the half space Vj C M"* corresponding to (Wlb related to the 
polytope 'PN,d^ Lemma[r|states that 



and 



VN,d c V, 



(42) 



(43) 



where the star * denotes here the restriction to spectra, i.e. 
ordered and normalized vectors. There are two possible re- 
lations between Vj (or Vj\^) and VN.d- They are illustrated 
in Figure |2] in form of a simplified 2— dimensional picture: 
There, we consider two half spaces Vi and V2 correspond- 




FIG. 2: Polytope "Pjv.d and two restricted generalized Pauli con- 
straints Di, D2 > with boundaries 5*1, S2 arising from two con- 
straints Di , D2 > belonging to a higher dimensional marginal 
settings A^['Hd']. 

ing to the "restricted" constraints D[ > and D'2 > with 
boundaries 5*1 and S2 and orientation indicated by stripes. 
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Such hyperplanes can either contain a facet of maximal (ex- 
ample Si) or lower dimension of Vn^ or they lie outside of 
Vn.cI (example S'2). The third case of a proper intersection 
is not possible due to Lemma |T| Every constraint D' with 
boundary S of its restriction D' lying outside of 'PN,d is a 
constraint, which is irrelevant for the pinning analysis since it 
has the form 



D'{X)=c + Di\"^) + 0{Xd+i), 



(44) 



where £)(A*'^) > is a constraint of the setting A^[Hd] with 
a boundary shown in Figure [2] as hyperplane ^2 and c > 
is some offset. Hence if the spectrum decays sufficiently fast, 
constraint D' is not saturated at all due to the offset c and 
thus irrelevant. Moreover, for every Pauli facet of Vpf^d corre- 
sponding to some constraint D > 0, Lemma [T] guarantees the 
existence of a constraint D' > in the larger setting whose 
projection D' coincides with D. We summarize these insights 
by stating 



Lemma 2. Given two marginal settings A^[Hd\ '^nd 
/\^\Hd'] with d < d' E H. Every generalized Pauli con- 
straint D' > of the setting A^[Hd'] relevant for the pinning 
analysis is given by a linear modification of some generalized 
Pauli constraint D > of the setting A^[Hd]> 



D'{X)=D{X''-)+0{Xd+i). 



Finally, we remark that for the important case d' infinite ef- 
fectively the same results holds but one has to deal with one 
subtlety. Since Vn.oo is described by infinitely many con- 
straints Lemma [T] guarantees for every constraint _D > of 
the setting A^[Hd], only the existence of a sequence of con- 
straints D'j > whose restrictions Dj > converge to the 
constraint D > 0. This means that condition ( |45| ) in Lemma 
|2] holds up to a small error e. 



D'^{X)^s + DiX'') + OiXd+i), 



(46) 



which can be made arbitrarily small by choosing appropriate 
constraints D'^. Hence, to minimize the technical effort we 
assume in our work that Lemma |2] holds in its original form 
also for the case d' infinite. 



stated, which \'^) in that case satisfies: 

D\^) = (koU + KiNi + ... KdNd) I*) = 0. (48) 

Applying this condition to the expansion of | ^t) in Slater de- 
terminants induced by the natural orbitals. 



(49) 



it implies Klyachko 's selection rule, which states that when- 
ever 



D\\) ^ 0, 



(50) 



the corresponding coefficient Ci vanishes. To show the 
strength of this selection rule we study states in the Borland- 
Dennis setting. The corresponding Hilbert space A^\Hq] has 
dimension (3) 20 and the generahzed Pauli constraints read 

iHTi 



Ai + Ag, A2 + A5, A3 + A4 < 1 (51) 

D^s) := 2 - (Ai + A2 + A4) > 0. (52) 

The normalization together with the non-negativity of the 
eigenvalues leads to 



1. 



(53) 



(45) Hence the constraints in ( 5 1 1 are always saturated and this 



ito 


(|48]i 






(Id 


-Ni 


- Ne) 1*) 





(Id 


-N2 


- N5) 1*) 





(Id 


- A^3 


- Ni) 1*) 


= 



(54) 

Klyachko's selection rule applied to ( [54| implies that every 
Slater determinant showing up in the ansatz (49i for j^*) is 
built up by natural orbitals with one index from each set 
{1,6}, {2, 5} and {3,4}. Those are the 8 states 1 1,2,3), 
|1,2,4), |1,3,5), |1,4,5), |2,3,6), |2,4,6), |3,5,6) and 
1 4, 5, 6). If the constraint p2| ) is also saturated the selection 
rule restricts this family of Slater determinants to the three 
states |1, 2, 3), |1, 4, 5) and |2, 4, 6) and in that case we find 



Ivl-a) =a|l,2,3)+/3|l,4,5)+7|2,4,f 



(55) 



C. Selection Rule.- 



D. Quasi-Pinning and modified Selection Rule. — 



In this section we state a selection rule which explains how 
the structure of the A^— fermion state j'l') e A^[Hd] simpli- 
fies if the spectrum of the corresponding 1— RDO is pinned 
to some Pauli facet of V^.d- Moreover, we apply it for the 
setting A^iHe]. 

Let's consider a state |^) with natural occupation numbers 
A = {Xi)f^-^ saturating some generalized Pauli constraint 



D{X) = Ko + Ki^i + • • • KdAd > 0. 



(47) 



In ifTSl . by introducing the creation and annihilation operator 
a\,,ak of a fermion in the natural orbital |fc) and the parti- 
cle number operators Nj. = a^a^, an important condition is 



In this section we show for the Borland-Dennis setting that 
any state j^*) e A'^['H6] whose natural occupation numbers 
are approximately saturating the corresponding generalized 



Pauli constraint ( 52 1 also fulfill approximately condition (48 1. 
We also quantify this relation. This result then guarantees that 
our Hartree-Fock extension will work for systems exposing 
strong pinning. 

As a warm-up and since we will need the result we first 
study a simpler question. It is a basic fact that the spectrum 
Asi = (1, . . . , 1, 0, . . . , 0) can only arise from a Slater deter- 
minant 1 5*) = |1, . . . , A^). Is this statement stable under small 
deviations, i.e. A w Agj I^I^) w |1, . . . , A^)? Yes, it is true 
according to 
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Lemma 3. Consider a state \^) £ A^lHd], let {\k))'l^^ be 
its natural orbitals and denote the projection operator onto 
the space spanned by |1, . . . , N) by Psi- Then, 



N 



where 



< N - {Xi + ... + Xn) =■■ 5. 



(56) 



(57) 



Proof. We expand the state j^*) in Slater determinants induced 
by natural orbitals (recall Section [A|, 



i*) = E 



We define the operator 



S = Nld - [a[ai 



(58) 



(59) 



Since all operators a|aj, i — 1, . . . ,d commute it is clear that 
S has the spectrum {0, 1, . . . , N} with eigenstates |i). The 
eigenvalue corresponding to |i) is the number of indices k Gi 
not belonging to the set {1, . . . , N}. We denote the set of 
indices leading to the eigenvalue k by Jk and find 

5 = - (Ai + . . . + Ad) 

= (*|Arid- (aIai + ... + 4ad) |*) 

= mm 

^ c;ci(j|^|i) 

ije,/oU...uJiv 

= E l^il'(il^li)- (60) 



ie,7oU...UJjv 

Since for i G Jfc, 



we find 



N N 

k=oieJk fe=iieJfe 



(61) 



(62) 



and alternatively also 



N 



fc=i ieJfc 

The normalization of X^ikiP 
{(1, 2, . . . , N)} contains only one element) 

^ '1-|C(1....,7V)|'<^. 



(63) 



1 yields (Jq 



< 



and thus 



N 



l-'5< |c(i,...,iv)l' <l--S. 



1 

N' 



(64) 

(65) 

□ 



Now, we come back to the original question. We first state 
the mathematical result and present the proof afterwards. 

Theorem 4. Given a state \'^) e A'^['H6] with natu- 
ral occupation numbers {Xk)\^i- Let P be the projec- 
tion operator onto the subspace spanned by the states 
|1,2,3), |1,4,5), |2,4, 6), which corr espo nds to exact pinning 
of D^^' = As + Ag — A4 > (recall (55 \). Then as long as 



i5 = 3 — Ai — A2 — A3 < 



(66) 



( which means nothing else but being not too far away from 
the spectrum Asi — (1, 1, 1, 0, 0, 0) of a single Slater determi- 
nant) we find 



1 



Xs 



^(6) < IIP^-II^ < 1 



1 



with 



XS 



1 + 25 
1-4(5' 



(67) 



(68) 



Proof In Section[c]we concluded that has the form 

I*) = a|l,2,3)+/?|l,2,4)+7|l,3,5) 
+ 5|2,3,6) +z^|l,4,5) +/x|2,4, 6) 
+ C|3,5,6)+C|4,5,6), (69) 

with natural orbitals {|fc)}f.^i. Since the corresponding 
1-RDO is diagonal w.r.t. 



{k\pi\l) = Ski Afc, 



we find (recall ( pTj l) 



A4 = W\' + W\' + \fA' + \C\' 

A5 = |7|2 + |^|2 + |^|2 + |^|2 
Ae = |J|2 + |^|2 + |^|2 + |^|2 



(70) 



(71) 
(72) 
(73) 



The goal is now to show that the coefficients /3, 7, S, £, and C 
are small, i.e. 



\P^\ 



L2 = |a| 

= l-(|,9p + |7p + l<5p + leP + ICP)(74) 



m 

|2 



is close to 1, whenever ( |52| ), which here reads 

d(^'^ = -m' + \^\' + \s\^+2\e + \c? 

is approximately saturated. First we observe 
1 



(75) 



pm 



< 1- 



< l-2(-|/3p + l7l' + l'5|' + 2|e|^ + |Cl') 



(76) 



which is the upper bound for ||Pvl/||2 ^ in Theorem|4] 
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To derive the lower bound note the essential difference in 
(|75 I and (76i, the sign of the term |/3p. To get rid of this 



we write |^ = -x + (1 + X > and estimate 

(1 + x)|/3p in terms of |7|2, \5\^, \^\^, \(\^. For this observe 



that ( 70 1 in particular implies 



(77) 



which leads by the triangle inequality, the identity (A + B 

C)2 <S{A^+B^+ C2) and \fi\^, ji/p, j^p, \C\^ <l-\a 
to 



1 



< 



< 



< 



a 

A i\i\W\ + \s\\f^\ + \^\\C\f 

{\l\'w\'+\s\'\^^\'+\em 



3(1 -|«| 



h\' + \5\' + ^{m' + \c\')](m 



Now, for all s, r > we find by using ( [78] l 

l/sp + M' + H' + ieP + icP 

< {l-rM' + \j\' + \S\' + {l + s){2\e + \Cn 



(1-r) 



1 



3r(l - lal 



(79) 



By choosing 



2|a|- 



4|a|2-3 
4(1 

4|a|2-3 



(80) 
(81) 



the last expression in (79 1 coincides with Z?*^^' up to a global 
factor X- Both parameters r, s are non-negative as long as 
|ap > |. Finally, this leads to 

\\pnh = ^-m' + h\' + \s\' + \e + \c\') 

> 1 - (r- 1)L>(6) 

= l-Xi-\a\-D^^\ (82) 



with 



= r - 1 = 



3~2\a\' 
4|a|2-3 

1 + 2(1-1^1^) 
l-4(l-|aP)- 



(83) 



Lemma [3] states jap > 1 — S and since x is monotonously 
increasing, Xi-|aP ^ Xi5' which finishes the proof. □ 
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